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Abstract. We establish metastability in the sense of Lebowitz and Penrose under practi- 
cal and simple hjqjothesis for (families of) Markov chains on finite configuration space in 
some asymptotic regime, including the case of configuration space size going to infinity. By 
comparing restricted ensemble and quasi-stationary measure, we study point-wise conver- 
gence velocity of Yaglom limits and prove asymptotic exponential exit law. We introduce 
' soft measures as interpolation between restricted ensemble and quasi-stationary measure 

I to prove an asymptotic exponential transition law on a generally different time scale. By 

using potential theoretic tools we prove a new general Poincare inequality and give sharp 
' estimates via two-sided variational principles on relaxation time as well as mean exit time 

and transition time. We also establish local thermalization on a shorter time scale and give 
mixing time asymptotics up to a constant factor through a two-sided variational principal. 
\^ ' All our asymptotics are given with explicit quantitative bounds on the corrective terms. 



1. Metastability after Lebowitz and Penrose 



1.1. Phenomenology and modelization. Lebowitz and Penrose characterized metastabZe 
thermodynamic states by the following properties ||3]| : 

(a) only one thermodynamic phase is present, 

(b) a system that starts in this state is likely to take a long time to get out, 
^ ', (c) once the system has gotten out, it is unlikely to return. 

^ ! We can think, for example, to freezing fog made of small droplets in which only one phase 

is present (liquid phase) that remains for a long time in such a state (until collision with 
ground or trees, forming then hard rime) and that once frozen will typically not return to 
ff-^ I liquid state before pressure or temperature have changed. 

O ■ To modelize such a state they considered in [ 3 ] a deterministic dynamics with equilib- 

rium measure fi. First, they associated with the only one phase of the metastable state a 
subset TZ of the configuration space, and described this metastable state by the restricted 
ensemble fin = fj,{-\TZ). Second, they proved that the escape rate from TZ of the system 
^ i started in fi-ji is maximal at time t = 0, and that this initial escape rate is very small. Last, 

^ ■ they used standard methods of equilibrium statistical mechanics to deal with (c) . As esti- 

- - - mate of the returning probability to the metastable state they used the fraction of member 

of the equilibrium ensemble that have configuration in TZ and they noted ([3], Section 8): 

This amounts to assuming that a system whose dynamical state has just 
left TZ is no more likely to return to it than one whose dynamical state was 
never anjwhere near TZ. The validity of this assumption, at least in the 
short run, is dubious, but at least it provides us with some indication of 
what to expect. 
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In this paper we want to give a different model for the same phenomenology that over- 
comes the last difficulty. We will work with stochastic processes rather than determinis- 
tic dynamics, but the Lebowitz-Penrose modelization will be our guideline. We will try 
to recover this phenomenology under simple and practical hypotheses only. Since the 
study of metastability has been considerably enriched after Lebowitz and Penrose work, 
we want also to incorporate in our modelization as much as possible of what was previ- 
ously achieved. We will then make a brief and partial review of these achievements. On 
this review will depend our goals and starting ideas but not our proofs since we want 
to make this paper as self-contained as possible. Our model and results are presented in 
Section m 

1.2. A partial review. Since Lebowitz and Penrose paper, an enormous work has been 
done to describe the metastability phenomenon. In particular Cassandro, Galves, Olivieri 
and Vares introduced the path-wise approach, which focused, in the context of stochastic 
processes, on time averages associated with an asymptotic exponential law [5]. This was 
further developed by the pioneering works of Neves and Schonmann who studied the typi- 
cal paths for stochastic Ising model on a given volume in the low temperature regime (lIZl], 
llSl]). This work was then extended to higher dimensions, infinite volume and fixed tem- 
perature regimes, locally conservative dynamics and probabilistic cellular automata r ill2ll . 
IfTm . IfTSll. 112311. m, IHOII). 

As developed in (24] , a crucial role was played by large deviation tools inherited from 
Wentzell and Freidlin in their reduction procedure from continuous stochastic processes 
to finite configuration space Markov chains with exponentially small transition rates [SJ. 
This is especially true for very low temperature regimes, but the same kind of reduction 
procedure allowed to deal in various cases with large volume rather than low temperature 
limits (see [,l4i1 for Curie-Weiss model under random magnetic field, see [241 for further 
examples). 

Then, using potential theoretic rather than large deviations tools, Bovier, Eckhoff, Gay- 
rard and Klein, developed a set of general techniques to compute sharp asymptotics of the 
expected value of asymptotic exponential laws associated with the metastability phenome- 
non, and revisited (after Q, IITOll ') the relation between generator spectrum and metasta- 
bility ([16], [19], [221). This allowed, for example, to go beyond logarithmic asymptotics 
for stochastic Ising models in the low temperature regime ( [,17i] . [|25,] ) and to prove the 
first rigorous results in the fully conservative case ( |[35lD . to deal with metastability for 
the random hopping time dynamics associated with the Random Energy Model ( IflSll '). to 
make a detailed analysis of Sinai's random walk spectrum ([28]), to extend the study of the 
disordered Curie-Weiss model to the case of continuous magnetic field distribution ( |[30ll . 
II36I1). 

We then reached to an essentially complete comprehension of the metastability phe- 
nomenon in at least two classes of models: very low temperature dynamics in finite fixed 
volumes and large volume or continuous configuration space dynamics that can be re- 
duced via Wentzell-Freidlin procedure to the previous case. Of course, specific and often 
non trivial computations have to be made for each specific model, but there exists a general 
approach to the problem that is developed in lf24|] and, as far as the potential theoretic part 
is concerned, llT6ll . |[T9l] . and |[22ll together with I138L that bridges between potential the- 
ory and typical path description by reinforcing and generalizing the results of 1[9]] (and it 
is worth to note that [38], after [i33li37 ], contemplates also the case of polynomially small 
rather than only exponentially small transition probabilities) . For both classes of models, 
like one-dimensional metastable systems as considered in |[28ll or II27I1 , a recurrence prop- 
erty to a very localized subset of the configuration space (single configurations identified 
to metastable states in the first case, small neighborhoods of the dynamics attractors in the 
second case) plays an important role. 
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Beyond these two classes of models there are many limit cases, special cases, and partial 
results. For example, in IfTSlI we are far from a finite fixed volume situation but single 
configurations can still be identified with metastable states and have still enough mass at 
equilibrium for potential theoretic or renewal techniques to work. This is not the case 
in [351 . where potential theoretic tools give only expected values of some hitting times 
when the system is started from some specific harmonic measures that are very different 
from what one would expect to be a "metastable state" (here, like in the sequel and fol- 
lowing Lebowitz and Penrose, we mean a whole measure when referring to a metastable 
"state" and not to a single configuration of the configuration space) . Any kind of exponen- 
tial law is presently also lacking in this case. The same difficulty is faced in [30], but it is 
overcome in [36] by mean of a specific coupling argument that gives point- wise estimates 
and opens the way to the exponential law. Finally, the beautiful paper by Schonmann 
and Shlosman IfTS]] achieves the tour de force of using essentially equilibrium statistical 
mechanics computations to deal with the dynamical problem of metastability. In this case 
also the exponential law is lacking as well as sharp estimates on the transition time, even 
though the simple formulation of such properties is not completely obvious in this fixed 
temperature and vanishing magnetic field regime. 

1.3. Starting ideas. In the present paper we want to elaborate some tools to describe 
the metastability phenomenon beyond the case of a dynamics with a recurrence property 
to a very localized subset of the configuration space. We will focus on exponential laws 
and sharp asymptotics of their expected values. We note that the exponential law itself 
suggests some kind of recurrence property. If it is not a recurrence property to a very 
localized subset, it has to be in some sense a recurrence property to a whole "spread 
measure". And this measure should coincide with our metastable state. Now, following 
Lebowitz and Penrose, if we associate the metastable state with some subset TZ of the 
configuration space X, then, considering property (b), we have at least two candidates to 
describe our metastable state: one is the restricted ensemble fin = /u(-|7^) the other is the 
quasi-stationary measure 

= lim G ■\T;^\n > t) (1) 

where X{t) is the configuration of the system at time t and Tx\n is the exit time of TZ (we 
will be more precise in the next section). The main advantage of fiji is that fiji is often 
an explicit measure one can compute with, while /^i^ is only implicitly defined. The main 
advantage of is that the exit law of TZ for the system started in is an exponential 
law. Our first results will then start with a comparison between fiji and We will give 
simple and practical hypotheses to ensure that they are close in some sense, then we will 
be able to prove some kind of recurrence property to fi*j^. All this will be done in the 
simplest possible setup: considering a Markov process on a finite configuration space in 
some asymptotic regime (including the possibility to send to infinity the cardinality of the 
configuration space) . 

In the present work we will essentially build on the ideas of four different papers: ^ 
for the formulation of the problem, [ 5 ] for the focus on exponential laws, fT6l for the 
introduction of potential theoretic techniques in the metastability field to get sharp esti- 
mates on some mean hitting times, and Miclo's work [34] where some concepts of local 
equilibrium, and "hitting times" of such equilibriums, are introduced. As far as this last 
paper is concerned, it will only work as a source of inspiration: we will not require a full 
spectrum knowledge, and we will not introduce any notion of a starting point depending 
local equilibrium. Finally, we note that the idea of considering quasi-stationary measures 
as metastable states was already contemplated in [20] . Even though some of our results 
echo some of [ 20 ] , we were not able to make any clear comparison, essentially because of 
the much more anal5^ical point of view of ||20]| and the many hypotheses introduced in the 
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results of Il20ll . We note that, on the one hand, Il20ll deals with a much more general setup 
than ours since the authors consider non-reversible Markov processes on a continuous 
configuration space, while we look at reversible Markov processes on finite configuration 
space. However, the reason why we assume reversibility is to be able to use potential the- 
oretic results to get sharp estimates on mean times via variational principles, a question 
that is not considered in Il20ll . 



2. Model and results 

2.1. Quasi-stationary measure and restricted ensemble. We consider a continuous time 
Markov process X on a finite set X with generator defined by 

cf{x) = Y,pi^^y)ifiy)-fi^)) (2) 

for X in X and f : X ^ M, and where p is such that ^yP{x, y) = 1. Since X is finite, any 
generator can be written like in dU up to time rescaling. We assume that X is irreducible 
and reversible with respect to some probability measure fi, we denote by (•, •) the scalar 
product in i'^ifj-), by || • || the associated 2-norm, by V the Dirichlet form defined by 

-Dif) = if, -cf) = \Y1 y) [/(^) - /(y)]' 

x,y£X 

where each conductance c(x, y) is equal to 

c(x,2/) = (4) 

and by 7 the spectral gap 

7 = mm — — -. (.5 J 

Var^(/)5^o Var^(/) 

For TZ a X ^Ne define in each x gTZ the escape probability 

and we denote by Xji the reflected process with generator given by 

yen 

for X in 7^ and / : 7?- — )■ M, and where, for all x, y in TZ, 

I \ f p{x,y) iix ^ y, 

Pn{x-,y) = { ) i /X -r (8) 
F/<.v I p{x,x) + e-jzix) ifx = y. 

We will only consider subsets TZ such that both Xn and X;^;\n are irreducible and we note 
that Xfi inherits from X the reversibility property with respect to the restricted ensemble 

fin = K-\n). (9) 

We identify ('^{fin) with the subset of f{fi) of functions f : X ^ R such that flxyii = 
and we denote by (•, \\ ■ \\ti, V-ji, c-jiix, y) and 77^ the associated scalar product, 2-norm, 
Dirichlet form, conductances for x, y inTZ and spectral gap. 

We denote by p^ the sub-Markovian kernel on TZ such that, for all x, y in TZ, 

Pn{x,y) =p{x,y). (10) 
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We know from [2] and the Perron-Frobenius theorem that there exists c/)^ > such that 
1 — (f>^is the spectral radius of and that there is a unique quasi-stationary measure /i^ 
such that /U^p^ = (1 — In addition we have, for all x,y inTZ and t > 0, 

limj^+oo Px{X{t) = y\Tx\n >t) = ^^(y), (11) 

P^,'^{Tx\n>t) = e-'f'V, (12) 

l^*nM = rn- (13) 

In Sections [2.31 and [6l these properties will be rederived in a slightly more general context. 

Our first result states that if 1/0^, the mean exit time for the system started in /x^, 
is large with respect 1/77^, the relaxation time of the reflected process, then the quasi- 
stationary measure /i^ is close to the restricted ensemble fxn. More precisely, defining, for 
all X in TZ, 

e*n = ^ (14) 

m 

hUx) = ^ (15) 

we prove the following. 
Proposition 2.1. Ife^ < I, then 

Var^^(/i^) = \\h*T^ - Infn < (16) 



1 — F* 



Proof. See Section [3TT1 □ 

Remark. When proving that goes to in some asymptotic regime (for example when 
the cardinality of the configuration space goes to infinity like in [5]] , when some parameter 
of the dynamics goes to like the temperature in [7] or when both happen like in Il32ll ) 
one has to give upper bounds on 0^ and lower bounds on 77^. As far are such upper 
bounds are concerned we will show that they can be easily obtained using test functions 
in some inequalities (see Theorem 12.111 below together with the following remark) . In 
addition, since one can often easily compute with ^-ji and en is often explicit, one can 
usually estimate 

(l^n = fJ-nien) (17) 
and use the following lemma that we prove in Section [3^21 

Lemma 2.2. cp^ < cpn. 

Lower bounds on 77^ can be more difficult to obtain. However we note, first, that rough 
lower bounds will often be sufficient to our ends, second, that the new Poincare inequality 
we will prove in this paper (Theorem 12.1 2D can be used to this purpose. 

As a consequence of this first result we can control the convergence rate of the Yaglom 
limit (fTTl) . We note that, by the reversibility of X with respect to fi, is a self-adjoint 
operator on i'^i^n) and has real eigenvalues. By the Perron-Frobenius theorem, this implies 
the existence of a spectral gap "j^ > equal to the difference between the first and the 
second largest eigenvalue of p^. 

Proposition 2.3. If < |, then 

Proof See Section [3T3l □ 



'<^l^}. as, 



METASTABILITY AND QUASI-STATIONARY MEASURES 



6 



Remark. Since, after the statical study made in Il39ll , we intend to apply our results to the 
dynamical study of the cavity algorithm introduced in Il26l] , for which finite effects volume 
are of first importance, we need to give asymptotics with quantitative control of corrective 
terms. This produces quite long formulas and to simplify the reading we put between 
graph parentheses any terms that goes to 1 in a suitable asymptotic regime. 

Then we define 

Cn = inmnnix)h^'^ix) = min/i^(x)/i^(x) (19) 
and, if < 1, for any 5 g]0, 1[, 

to get point-wise mixing estimates for Yaglom limits. 

Theorem 2.4 (Mixing towards quasi-stationary measure). If e"^ < 1/3^ then for all x,y gTZ 

and 5 g]0, 1[, 

F,{Xit)=y\Tx\n>t) 



<5 as soon as t > Tl^^ (21) 



Proof See Section [3!4l □ 

Remark. In words this says that either the system leaves IZ before time T^*^ or it is de- 
scribed after that time by /^^ in the strongest possible sense. This theorem is useful only 
if one can provide upper bounds on T^-j^. We already estimated 7^ in Proposition 12.31 
and we note that, as far as C,^ is concerned it only appears in the formula through its 
logarithm. Very crude estimates will then often be sufficient. One has for example the 
following lemma. 

Lemma 2.5. Ifp{x,x) > 0/or all x eTZ, then 

C^>mmf,*T^\x)>e-''^^''^, (22) 

where 

Ati = max{- lnp7^(x, y) : pn{x, y) > , \/x,y e TZ} 
Dti = min{k > : Pj^ix, y) > , \/x,y e TZ} 

Proof See Appendix |Al □ 

We will make a special choice for the parameter 5 in (l20l) . We define 

= (n e-jiix) (23) 

and set 

= n*^,n- (24) 
The next lemma is necessary for the following results to make sense. 

Lemma 2.6. For e\ < I, it holds 



Proof See Section HTTl □ 

Remark. We will often refer in the sequel to the regime ^ 1. The lemma implies that 

in this regime, 6^ and (/)^ also go to 0. In addition, it provides a sufficient and practical 
condition for being in such a regime: since minx^a^n e-jiix) is usually an easy to estimate 
quantity, one has a lower bound on 6^ as soon as one can estimate from below, and, 
since on the right hand side 5^ appears only through its logarithm, rough bounds will often 
be enough to show that goes fast enough to 0. 
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We close this section with a first asymptotic exponential law. 

Theorem 2.7 (Asymptotic exit law). For any probability measure v on TZ, define 7r7^(z/) 
F.(TA\7e < T^)- lfe*T^ < 1/3, then 
(i) 

/ E,[r;^n]<T^+'-^{l + S*n}<i^{l + 4>*nTn + ^*n} 



(ii) for all t > (t)^T^, 

r F,{t;^t^ > ^) < (1 _ nn{i^))e~' {eW{l + <5^)} 
1 Kirx^n > ^) > (1 - 7Tn{u))e-' {6^(1 _ 5*^)} 



Proof See Section [4^21 □ 

Remark. The theorem gives more than asymptotic exponential exit law. The first part 
shows in particular that in the regime 0^T-^ ^ 1, the quasi-stationary measure maximizes 
asymptotically the mean exit time on the set of all possible starting measures. The second 
part shows that, provided it'ji{u) converges to some limit, the normalized mean exit time 
4>Ti^?(Sji converges in law to a convex combination between a Dirac mass in and an 
exponential law with mean 1 . 

As an example of application we can consider the case of the restricted ensemble. 

Lemma 2.8. It holds 

TTnM = F,.^(T;t\7e < T^) < <PnT^ ■ (26) 
Proof See Section RTSl □ 
Corollary 2.9. If < 1/3, then 

^ {(1 - <PnTm + ^nT^ " ^ < ^^^n[rx\Ti] < ^ {1 + ^nT^ + ^*n} (27) 
and, for all t > <p^T^ , 

e-* {(1 - ^nTm - S*n)e'''^''^} < ^,Arm > J^) < e'' {{1 + 5^)e'^-^-} . (28) 

This shows asymptotic exponential exit law in the regime (pnT-^ ^ 1 for the system 
started in the restricted ensemble. 

2.2. (k, A)-capacities, mean exit times and a new Poincare inequality. In this section 
we introduce a new object which extends the notion of capacity between sets. For any 
K, A > and A,B c X, we first extend the electrical network {X, c), with c{x, y) = 
IJ,{x)p{x , y) = fi{y)p{y,x) for all distinct x,y G X, into a larger electrical network {X,c) 
by attaching a dangling edge (a, a) with conductance «;/u(a) to each a e A and a dangling 
edge {b,b) with conductance \fi{b) to each b e B. More precisely, we add \A\ + \B\ nodes 
and edges to the network by setting 

X = X\j{a : aeA}\j{b : b £ B} 

and, for all distinct x,y G X we define 

c{x, y) if (i, y) = {x,y) e X X X 

Kn{a) if {x, y) = (a, a) for some a £ A , , 

\filb) if {x, y) = {b, b) for some beB ' 
otherwise 



c{x,y) 
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Definition 2.10. The (k, A)-capacity, C^{A, B), is defined as the capacity between the sets 
A and B in the electrical network {X,c), and then is given by 



min \\ c(x,y)[/(x)-/(y)]2;/| =1,/| =oi 
min {V{f) + K ;u(a)[/(a) - l]^ + A J] ^(6)[/(6) - 0] 



(30) 



= f^^R ^^^^^ + '^/«(^)IE,M [(/U - 1)1 + - 0)'] } • 

Remarks. 

i) Since all the points of A and B are at potential 1 and respectively in formula (|30]) . 
they are electrically equivalent and we could have defined the {k, A)-capacity be- 
tween A and B by adding just two nodes to the electrical network {X, c). However, 
our definition with dangling edges will be more useful in the sequel. 

ii) A (k, A)-capacity is in some sense easy to estimate since it satisfies a two-sided 
variational principle. On one hand, by definition, it is the infimum of some func- 
tional, and any test function will provide an upper bound. On the other hand 
it is the supremum of another functional on flows from A to B, which are anti- 
symmetric function of oriented edges with null divergence in X, i.e., on functions 
4) : X X X R such that for all x e X \ {AU B), div^V^ = Exe;? = 0- 
This is Thomson principle that goes back to [1], Chapter 1, Appendix A (see also 
lecture notes llSTll for a more modern presentation or textbook I113II for the proof 
of an almost equivalent result) . Letting 



be the energy dissipated by the flow -tp in the network {X, c), and 4)i{A, B) the set 
of unitary flows from A to B, that is, the set of flows -ij) from Ato B such that 



^ diva^ 



v;(a, x) = 1 = - ^ div-,v; = -Y.Y. ^(^' ^) ' <^3i^ 



b&B 



or, if P(^) < +00, such that 



(32) 



b&B 



we have 

C^iAB) 



max _ T>{'ip) 
i>&4>i{A,B) 

max i V{%p) + y — h > , 



(33) 



max^^<|p(V) + ^E,,^ 



where (t>i{A, B) is the set of unitary flows il) from Ato B and 

'4){x,y) 



/'(g) F 



(^) 



c{x,y) ' 



Then, any test flow provides a lower bound on C^{A, B). 
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iii) We know (IfTSH. I13T1D that the infimum and supremum in (l30l) and (l33l) . are 
reahzed, respectively, by the equiHbrium potential = P(.)(Ga(o"k) < Gs(crA)), 
where Ga and G b are the right continuous inverses of the local times in A and B, 
while (Tk and cja are independent exponential times with rates k and A, and by its 
associated normalized current 

- OT) ■■{^,y)^X^X^ ^^^iy^i^:,) _ v^{y)) . (34) 

We will say more on such quantities in the next section. 

iv) The previous definitions and observations extend to the case when k and A are 
equal to +00. In that case we identify A to A in the extended network if k = +00, 
or i5 to i3 if A = +oo,and we drop the infinite upper or lower index in the notation, 
so that, for example, C,^{A,B) = C'i^{A,B). However, when k and A are both 
equal to +00, to avoid any ambiguity we need to require that AD B = (/}. In that 
case the notation becomes C{A, B) = C'^{A, B) and we recover indeed the usual 
notion of capacity. 

We then get sharp asymptotics on mean exit times for the system started in the quasi- 
stationary measure. 

Theorem 2.11 (Mean exit time estimates). If < 1/3, then, for all k > 0, 

— jKtT) — y-<^n-^]^(Pn^ — jKji) — + </)^i^ + + — | ^.ibj 
Proof See Section [5l □ 
Remarks. 

i) In the regime + </>^T^ ^ 1 one can choose k such that cl)^<^ n <^ and get 
matching bounds. 

ii) Both bounds are in some sense easy to estimate since capacities satisfies a two- 
sided variational principal. Moreover, compared with the formula for mean exit 
time provided by potential theoretic techniques (see, e.g., [16]), the above in- 
equalities require no residual mean potential estimates. (Such estimates, as well 
as some harmonic measures will only play a role in the proof oi the theorem.) 

iii) Resolving the upper bound in cj)^ (0^ appears also in the right hand side of the 
inequality) one can get good bounds on cj)^ and prove that it decreases fast enough 
to guarantee 4>^T^ <^ 1. 

Our (k, A)-capacities provide also spectral gap estimates and a new general Poincare 
inequality. For k, A > and A, B c X we set 

Theorem 2.12 (Relaxation time estimates). For all /t, A > and any TZ C X such that X-ji 
and Xxsji are both irreducible Markov processes, 

Kfi(ll) A^^(AV^) / 

7 - </.^(7^,A^7^) V V ' 7A^7^ 

Proof See Section [5j □ 
Remarks. 

i) Without loss of generality, we can assume fi{TZ) < ^l{X\'R) so that, by (1361) . 
^re < 2C^{TZ,X\TZ)/fi{TZ). Then, as a consequence of the previous theorem and 
of the monotonicity in k and A of {k, A) -capacities, we get matching bounds over 
1/7 in the regime + £*p^ + (l>nhx<ji ^ 1- O^ie can indeed choose k such that 




9} 



(37) 
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(j)^ <^ K <^ jTi, as for Theorem [27TT] (Remark i)), and A such that 0^, ^ A ^ 
'jX'^R.- In addition and hke previously, all the relevant quantities can be estimated 
by two-sided variational principles. 

ii) The lower bound is a generalization of the classical isoperimetrical estimate that is 
recovered for k = A = +oo. 

iii) The upper bound is a new Poincare inequality. This inequality, or an easy to derive 
version when one divides the configuration space into more than two subsets, 
echoes Poincare inequalities given in [|21.] . We are not able to compare in full 
generality our result with that of I1 2T11 but we note that because of the presence 
of some global parameter called 7 in [21] one gets generally in our metastable 
situation an extra /actor 1/ mm{^'ji,^xxiz) ^PPlyii^g results of [|21|] . 

2.3. Soft measures, local thermalization, transition and mixing times. We address 
now the difficulty raised by Lebowitz and Penrose. Whatever the measure we choose 
to describe our metastable state, restricted ensemble or quasi-stationary measure, it is 
associated with some subset TZ of the configuration space. Then there is an ambiguity 
when one looks at property (b): what is "getting out" of the metastable state? One is 
tempted to say that it corresponds in our model to the exit of TZ. But doing so we are very 
unlikely to modelize in any satisfactory way property (c) : we can expect that "on the edge", 
when the system just exited TZ, it has probabilities of the same order to "proceed forward" 
and thermalize in X\TZ and to go "backward" and thermalize in TZ. Then we would like 
to define what would be a "true escape" from TZ. Theorem 12.41 suggests an answer in the 
regime (p*x\jiT^^ <C 1. We could define the true escape has the first excursion of length 
Tj^^ inside X\TZ. Since time randomization is almost always a good idea, we are led to 
the following definitions. 
For any A c X we call 

LA{t) = [ lA{X{s))ds (38) 
Jo 

the local time associated with A up to time t and we denote by Ga the right-continuous 
inverse of La ■ 

GA{t) = inf{s > : Lr{s) > t}. (39) 

For ax an exponential time with mean 1/ A that is independent from X we define for all x 
and y in 7^ 

p*T^^xix, y) = (X(r+) = y, L;t^7^(4) < ^a) (40) 

with the return time in TZ after the first clock ring r, i.e., = t + t-ji o 0^ with 9 
the usual shift operator. (One can see the process X as a process updated according to its 
discrete version with transition probability matrix p at each ring of a Poissonian clock with 
intensity 1.) We also define, for all x in TZ, 

en,x{x) = Px■(L;t^7^(4) > ^ta) = 1 - V P*n,x{x, y) (41) 



^/67^ 



and for all x and y inTZ 



The Markov process Xji^x on TZ with generator defined by 



Cn,xf{x) = PnAfiy) - fix)) (43) 
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is reversible with respect to ix-ji and has spectral gap 

T^n,xif) .... 
77^,A = ^^ mm — (44) 

where 

VnAf) = \Y.''nAx,y){f{x)-f{y)f (45) 

x,y 

with 

cn,\{x, y) = fJ'n{x)pn,x{x, y) = Pn,x{y, a;)/i7^(y)• (46) 

In addition we define 

TX\Ti,x = Lii{Gx\n{(^x))- (47) 
We know by the Perron-Frobenius theorem that the spectral radius of ^ is a simple 
positive eigenvalue that is smaller than or equal to 1 and has left and right eigenvectors 
with positive coordinates. We call it 1 — and denote by fi^ the unique associated left 
eigenvector that is also a probability measure on TZ. 

Lemma 2.13. It holds 

ii) P^^^ (r^^,;, >t) = e-*^^.^ , Vt > ; 

iii) lim P^(X o Gnit) = y I r;t\7e,A > t) = /"^,a , Vx, y G 7^. 

Proof. See Section [6?T1 □ 

We say that fi^^ is a quasi-stationary measure associated with a soften barrier, or a 
soften quasi-stationary measure, or, more simply, a soft measure. Indeed, ^ is the limit- 
ing distribution of the process conditioned to survival when it is killed at rate A outside 7?-. 
So, the hardest quasi-stationary measure associated with TZ, corresponding to A = +oo, 
is the quasi-stationary measure /x^, while the softest measure, corresponding to A = 0, is 
the restricted ensemble fin (0^ o ~ ^ o equilibrium measure associated with 

PnQ~ PT^fi' which is reversible with respect to fi-ji). More precisely we have the following. 

Lemma 2.14. The function A G [0, +oo] ^-^ ^tja ^ ^^(A^7^)* a continuous interpolation 
between the restricted ensemble fiji and the quasi-stationary distribution ^u^. In particular, 
for any Aq € [0, +oo] and y elZ, we have 

lim f^n,x{y) = A^7^,Ao(y) (48) 

A-s>Ao 

and for all x eTZ it holds the commutative limite property 

lim lim F^{XoGTi{t) = y \ Txyn,x > t) = lim lim F^{XoGTi{t) = y \ tx^^^x > t) • (49) 

A— 5>Ao t^oo t— >oo A— >Ao 



Proof See Section [6^ □ 

Analogously to what was done in the case A = +oo we set ^tix ~ '^n xh^^,^' ^nx ~ 
l^Tix/l^Ti and we call 7^ ^ the gap between the largest and the second eigenvalue of ^ 
(since p^ is self-adjoint with respect to (•, it has only real eigenvalues). We also define 

(l>n,x = fJ'n{eii,x)- 

Proposition 2.15. The parameters "fn^x, (t>Tix, £7^^ '^^^ (t'n,x depend continuously on A. In 
addition, when A decreases to 0, so do cf)^ ^, and (/i-ji^x, while ^ii,x increases. 

Proof See Section [6l3l □ 

The proofs of Sections [3] carry over this more general setup, and we get the following 
theorem. 
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Theorem 2.16 (Mixing towards soft measures). For all X > 0, cp'^x - 't'li' ^t^,>^ — '^^^ 



nx — ^Ti' Proposition \2. 1 1 Proposition \2.3\ Theorem \2.4\ and Lemma \275\ hold with an extra 
index A and writing X o Gjz instead of X. 

Remark. By continuity and monotonicity, the hypothesis < 1 and e"^^ < 1/3 are 
always satisfied for A small enough. 

We are now ready to deal with local thermalization. We define inductively, for k, A > 0, 
the stopping times for i > 0: 



0, 



(50) 
(51) 
(52) 



To = 

Then for 5 € (0, 1) we call iq the smallest i > 1 such that one of the two following condi- 
tions holds, 

i) X{Ti) G n and Lnin) - Lnin-i) > T;^ ^, (53) 



It 



X{Ti) TZ and LA^7^(7-^) - Lx\Tiin-i) > 



(54) 



and we set ts = . 



Theorem 2.17 (Local thermalization). For any 6 g]0, 1[ and any probability measure u on 
X, ife'^ x < 1/3 and e*^ ^ < 1/3, then Vx G X, 



max 



MxiTi) = x\xiTs)£n) 



F.iXin) = X I X{ts) 7^) 



<5. (55) 



(56) 



□ 



Moreover if i = max {e!^'^^-'^-^ - 1, e^^**'^'^'^ - 1^ < 1, it holds 

P.(T-5>t(i + i))<e-{^} 

Proof See Section [6l4l 

Remark. For k and A small enough, we have e^^^ < 1/3 and ^ ^/"^- Then, when 

K and A decrease to 0, we have non increasing upper bounds on T^^ ^ and T^p^ ^. As a 
consequence, the condition ^ < 1 will always be satisfied for k and A small enough and the 
theorem says that starting from any configuration the system is close to a random mixture 
of two states (/^^ and /i^^^ ^, close to fiji and fix^ji respectively) after a time of order 

As previously me make a special choice for the parameter 6 in previous theorems. We 



set 



^^,a = C^,aX mill e7^(x) 



^7^,A — -'5^^,7^,A• 



(57) 



Note that here we use min^^gg T^ 67^(2;) and not min^^g^ T^ eTi^x{x) for the sake of concrete- 
ness and practical estimates : there is no simple way to give lower bounds on eTi^x{x) while 
min^^gQ T^ eTi{x) is usually an easy to control quantity. In doing so we have no straight 
equivalent of Lemma [Z6l We only have the following result. 



Lemma 2.18. Ife^^^ < f^^" 



\TZ \ < 2e?j i ( In T« 



!^7^,A ^ '/^7^,A-'7^,A ^ ^^7^,A y-^j. 
Proof The proof runs similarly to that of Lemma [Z6l and we omit it, 



1 + 



(58) 

□ 
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Now the proofs Section [4] carry over this more general setup and we get asymptotic 
exponential laws for the transition time tx^^\ that is the time spent in IZ before "truly 
escaping" from IZ as indicated by Theorem l2.21l below. 

Theorem 2.19 (Asymptotic transition law). For all A > 0, (f)ii,x < (pn and Theorem \2. 71 
together with its Corollary hold with an extra index A. 

We can also give sharp estimates on the mean transition time and asymptotics of the 
mixing time. 

Theorem 2.20 (Mean transition time estimates). If < ^Z^' then, for all k, A > and 

T> ^ + ^,and setting c/)^ = cl)^{n,X\n), it holds 



> 



Ci(Tl,XKR.) 



I ("^"^)-^(7)/("^"^) a4 f 1 - max f ^ 



(59) 



where 



«2 = e" 
«3 = e-^^.^^^.^ - ^ 

04 = 03(1 + ^*Ti,xTn,x - ^n,x + 



A 



Proof See Section 1631 □ 

Remark. In the regime + + (pn/^xxiz ^ ^ assuming fi{TZ) < ^{X\TZ) one can 
choose K and A in such a way that 0^ ^ k ^ 77^ and 0^ , (^^7^ ^ A ^ 7a\r> choose T 
such that 

- + i < T < (60) 

AC A (p*^^^ 

and then get matching bounds provided 0^ ^^T^ a ^ ^- Oi^*^^ again, all the relevant quan- 
tities can be estimated via a two-sided variational principle. 

Theorem 2.21 (Mixing time asymptotics). ^/^^t^q'^tj a < '^1^' define T = Gx<ji{'^x) 
and, for any x in X we write = "^xi^iT) = •)• Then, 

W^x - ^^?^(^^\\TV < ^T*^ Q + -5*^Q, (61) 

Px - hWtv < /^W + ^T*^Q + (62) 

In addition, if 

V = + XT*^^ , + i (5^^ + 5^,, + % + '/'^,a7^^,a) < 1, (63) 



then, with 



xex 

we have 



Tmix = maxinf <; t > : Wx{X{t) = ■) - /i||iv < \ (64) 



Tmix — /I \ i* • (65) 

(1 - V)rn,x 

Proof See Section [6^ □ 
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Remark. The theorem makes sense when fi{TZ) < ij.{X\TZ) and in the regime + + 
'^n'^Xsjio ^ ^- '-^^ then choose A such that (p*^ ^ A ^ ^A\7^o' ^^^> P^vided 

'^n A ^ ^' ^^'^ parameter r] can be as close to ^jl{1Z) < 1/2 as we want. Then, since the 
spectral gap goes like cjy^ ■^/fi{X\TZ) and ^{X\TZ) > 1/2 the theorem provides the correct 
order for the mixing time. 

Let us finally summarize our results. To have a mathematical model of the metasta- 
bility phenomenon described by properties (a) -(c), we describe metastable states by soft 
measures associated with a subset 7^ of a finite configuration space X where runs a re- 
versible Markov chain with respect to some probability measure in the regime £ti+£*x^ + 
(\>\hx<jt ^ 1 and such that /i(7^) < fi{X\TZ). In this regime all soft measures are close 
to the restricted ensemble (Theorem 12. 16D . If we choose k and A such that (j)^ <. k j-ji 
and (j)^, (f)*^ X 1x\r. then we can show 

i) local thermalization towards the soft measure /i?^ a or nx\n,K starting from any 
configuration in X and on a short time scale ^ + j (Theorem 12. 17D . 

ii) exponential asymptotic transition time on a long time scale ~ c'^inx^) ("the- 



orems [ZT9l and [2^20]) , 



iii) return time to metastable state on a still longer time scale -r^ — ~ „f/iy'^!p, (The- 
orem[Z20] applied to A'\7e in place of 7^). 

In addition relaxation and mixing times are of the same order as the mean transition time 
(Theorems 12 . 1 21 and 12 .2 ID - in particular the relaxation time has the same exact asymptotic 
up to a factor fi{X\TZ) - while exit times are on long, but generally shorter, time scale 
(Theorem l2.11D . And we note once again, that all relevant quantities can be estimated via 
a two-sided variational principles. 

3. Analysis in I'^inn) 

3.1. Proof of Proposition l2. ll We recall that the dynamics on TZ with reflecting barrier on 
d-TZ, is reversible w.r.t. fin with spectral gap 77^. In particular, for any function / G i^ifJ^n), 
it holds the Poincare inequality Var^^(/) < -^V-j^^f), where V'ji{f) is the Dirichlet form 
of / given by 

VnU) = {f,-^nf)^ = Yl Mx)f{x){6,{y)-pnix,y))fiy). (66) 

Applying the Poincare inequality to /i^, we get 

Var^^(/i^) < —Vn{h\) = — IJ''R.{x)h'^{x) {6^{y) -pnix^y)) h\{y) 



= — [ fJ-nih^^) - ^ fJ'n{x)pn{x,y)h^{y) 
= — I A^7^(/^^^) - V ii\{x){j)(x,y)^b^{y)e'ji{x))h\{y) 



(67) 



' x,y&l 

where in the first two lines we use that, by definition of 

m{x)h\{x) = ^\[x) and i^^{h^) = i^nih^^ 
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and then we exploit the definition of pn and p^. From the last line, using that is a left 
eigenvector of p"^ with eigenvalue (1 — 0^), we get 

Var^,(/i^) < ^f^nih*n') = ^Var^^(/i^) + ^ . (68) 
77^ 77^ in 

Finally, rearranging the terms in the above inequality and from the hypothesis e^j = ^ < 
1, we obtain the required upper bound. 



3.2. Proof of Lemma [2121 Let us denote by £^ the sub-Markovian generator associated 
to the kernel p^. For any function / G i'^ifj.fi): this is defined as 

{cy){x) = -fix) + J2 Pni^^ y)fiy) ■ (69) 

We have 

07^ = mm ——^ . (70) 

Taking / = 1?^ as test function, we get 

(pn^Y^ f^nix) 1 - X] y) 
xen \ yen 

x&n \ yen 

= X] l^'R-{x)e'jz{x) = (1)11 , 
x&n 

which concludes the proof. 

3.3. Proof of Proposition [221 The second smallest eigenvalue of the sub-Markovian gen- 
erator cj)"^ + 7^, satisfies the variational formula 

<A^ + 7^ = min{^^^^-%^^ :/^0,(/,/i^)7e = o| 

I {l-,J)n J (72) 



(71) 



mm 



{ (/, -cy)n : (/, }A)n = , (/, /)7e = l} 



Let / be a function on IZ that realizes the minimum in the above definition, with (/, f)'ji 
1. Since (/, h^-ji = 0, we have 

(/, h\ - ln)n = -if, ln)n = -f^nif) 



f^Uf)<\\f\\n-\\h*n-^nfn<Y-^- ^^^^ 



and then, by the Cauchy- Schwartz inequality together with Proposition 12. 1[ 

_* 

-n 

Now, writing the orthogonal decomposition / = fj,n{f) + 9, with finig) = 0) we have 

l = \\ffn = ^^Uf) + \\9fn 

and thus, from ([73l) . 

P* 1 _ Op* 

Using 5 as a test function in 



^ ^n ^ ^n 



77e = min|^^ :/^^0,^7^(/i) = o| , (74) 
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we get 



77^ < ^-^ma) = ^^VnU) . (75) 



We will conclude the proof using the following result. 
Lemma 3.1. For all f e i'^ifJ-n), it holds 



Vn{f)<^^={f,-cy)n. (76) 



Proof of Lemma WA\ . For all x,y eTZ with x ^ y, pn{x, y) = p{x, y). Then we have 

-^n{f) = \ E Mx)PTi{x,y)[f{x)-f{yf 

I fJ-n{x)p{x,y)[f{x) - f{y)f , 



(77) 



2 



Vnif) <\Y. Mx)p{x,y) [fix) - f{y)f < ^ , (78) 



since only the terms in x / y matter in this sum. Thus 

^mx)p{x, y) [J yxj - J yyji ^ 

and this provides the stated upper bound. 

To prove the equality, we recall that the space i'^ifi-ji) is identified with the subset of 
functions / G ^^(/u) with /|^^ = 0. Since, for all x,y gTZ, it holds that ijlti{x) = ^{x)/ ^{TZ) 
and p^{x, y) = p{x, y), we have 

^)=J^)I1 Kx)fix){5Ay)-p{x,y))f{y) 

= t^nix)f{x)i5,iy)-pUx,y))fiy) ^^9) 

= {f,-cy)n, 

that concludes the proof of the lemma. □ 

Turning back to the proof of the proposition, and combining inequality (|75]) with Lemma 
13.11 we get 

m < ^^{f, -^nf)n = ^^^{<t^n + ik) , (80) 

since / was chosen in order to have (/, —C*^f)'ji = (I>ti + 7^- Setting cj)"^ = e^jn and 
rearranging the terms in the last inequality, we get 



m s - — — In 



that under the hypothesis < 1/3 implies 

1 {\-eX. 



< 



n 
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3.4. Proof of Theorem 12. 4[ The proof is based on a classical trick to control mixing times 
with relaxation times. For any probability measure v on TZ, any f : TZ ^ M. such that 
f^n(f) 7^ ^ ^^'^ ^5 * ^ 0' oils '^^^ check that 

E,[/(X(s + t))l|,^^>,+i|] - fiUf)^u{rx\n >s + t) 

Indeed, one can rewrite the right hand side of the above equality as the sum of four terms, 
two of which coincide with the two terms in the left hand side by the Markov property, 
while the other two terms cancel using the quasi-stationarity property, i.e. 



(81) 



E^l, [f\rx,n>t]=f^nif)- 

As a consequence one gets 

|E,[/(X(s + t))l|,^^>,|] -/.|e(/)P,(r;t\7e >s + t)\ 



yen 



X (E,[/(X(t))l|,^^>,}] -P,(r^^ > t)t,Uf) 
and, by the Cauchy-Schwartz inequality 

|E,[/(X(s + - M^(/)P.(r;t^7^ >s + t)\ 



< 



IT 



TZ 



n 



We now estimate these two factors. Noting that 

P,(X(.) = . , r;t^7^ > = i^e<(-) and ^^.)[f{X{t))l{r^^>t}\ = e*^-/(-) 
and diagonalizing the self-adjoint operator £^ in orthonormal basis, one gets 



/i7^(■) 



< 



n 



with 0^ G [0,vr/2[ such that cose,y = {^,h\) = v{h*^) , and 



E(.)[/(X(t))l|,^^>,}] - ||/|Up^cos0; 



< 



7^ 



2 sin2e.e-2*('^5^+7^) 

TZ J 



with 0/ G [0, tt] \ {7r/2} such that 
Moreover, since 



lPl'(TA\7^ > s) = ^in 
from inequality (|84l) we get 



I _u_ II COS 6*,^ 



< (l7^, 

< II — 





h* 






Pl'(T-A^7^>s) 


II II 






^sin Oi^e'" 


<t>n+ln) . 



(82) 



(83) 



(84) 



(85) 



'7^ 



(86) 
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Using inequalities JMl) and (|86l) . we finally get 

^ "^vK Tx\TZ > S)hTi{-)\\^ 

^ W>..(Xls\ = . T„ 



1P,(X(S) = -,TX^TZ>S) 



COS ( 



+ 



< II— 

— II \\n \ 



1- II^Tilk) sin6',, e 



he" 



7e 



1/ II cosdi^ -<f>i,s\ 1* 



n 



n 



which provides an estimate of the first factor in (l83l) . 
To what concerns the second factor, noting that 

^{■){rx\n >t)= E(.)[l7^(X(^))l|,^^>,}] 

and that, from the definition of cos 0/ applied to / = In, 



cos t>l^ 



and 



sm 



\\h* IP 



we get 



2 ^ Var^^^(fe;,) ^_2t(^^+^*) 
II^'-rIItj 



and from inequalities (l85l) and (l88l) . 



(87) 



(88) 



n 



< 
+ 



(rvvi? > t)ll/IL \\h%, |L cos 



cos^^f e"'^'^' 



7^ 



cosOte 



< ll/lksin^/e- 



l/l^ 



7z 1 1 ' 1 1 7^ cos t// e 



7^ (89) 



which provides an estimate of the second factor in ([831) . 
Inserting dSTl) and ([891) in Eq. dSS]), we then obtain 

\E,[f{Xis + t))l{.^^>t}] - f^WWuiTx\n >s + t)\ 
<\\^\U^ + \\h*n\\n)sm9, 



To conclude our proof we will make two more steps. First notice that from (l86l) one also 
get that, for any t > 0, 

P.(TA^7^ >t)> II^IU (p|te-*<^^ - sin0,e-*(<^^+^^) 



In particular, as soon as the following condition is verified 

sin0^e'*('^^+^^) < (J^^e--^^* ^ ||/i^|L tan 0^ e""^^* < 5, 
li tjIItc 

it holds 

IP.(TA^7^ > > (1 - 



(91) 
(92) 
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Now, dividing both terms of ( [901 ) by iJ'^{f)^v{Tx^ > s + t), we reach an inequality that 
controls the Yaglom limit and that, provided condition ( [9T1 ) and from the last inequality, 
reads as 



IE4/(^W)kA^7^M] 



< 



1-5 



(93) 



As a final step we apply this inequality to v = 5x and f = 6y. For this choice of u and /, 
and by definition of 9i, and Of, one has 



tan^i, < 



\h: 



and tan 6*/ < 



\h: 



Thus, from ( [93l ), we obtain that under condition ( [91] ) 



cos ( 



P^((X(t)=y)|r;t^7^>t) 



< e 



(1 + I|fc|^ll7z)ll^;z! 



+ 



/Var^^(/i* ) 



(94) 



+ 



l~e* 



where in the second line we used that > 1, the estimate given in Proposition 12. 1[ 

and we introduced the quantity defined in ( [T9l ) . 

The right hand side of the last inequality is smaller than 6 as soon as 



t > 



In 



5(1-5)C 



^ + ln 



1 + 



1 + 



(95) 



which also implies ( [91] ) . 

Finally, from the hypothesis < 1/3, from the concavity of the logarithm and of the 
square root function, and using that < 1, then 5(1 - S) < l/A and ln8 > 1 + 5/(4\/2), 
after some computation one obtains that the condition ( [95] ) is implied by the stronger 
condition 



t > 



In 



1 + 



(96) 



which concludes the proof. 



4. Around the exponential law 



4. 1 . Proof of Lemma 12.61 The proof uses the super-harmonicity of on TZ, where 
is thought as a function in ^^(/u) with h^{x) = for x / 7?-. To prove that is a super- 
harmonic function, notice that, for allx,z eTZ, 

-hnix) + ^p{x,y)h\{y) 



{Ch\){x) 



y&X 



-h*n{x) + lim 



p{y,x) ^,{X{t)=y,Tx\Ti>t) 



(97) 



i-5.oo ^-^ Unix) 



>t) 



-hnix) + lim V 



,{X{t + T) =X, TX\n >t + T) 



where t + r is the time of the first clock ring after t. Then 

F,{X{t + T)=x\TX\n>t + T) _ 



{Ch*T^){x) < -h^x) + }im J2 



l^n{x) 



hUx) + hUx) = 0. 
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This implies that /i^ reaches is minimum in TZ on d-TZ. As a consequence, being /i7^(/i^) = 
1, there exists xq € d^TZ such that h^{xo) < 1. Setting e^^ = min^g^-Te ^nix) and recalhng 
that = mm^fz-ji we then have 

= e7^ < e7^A^^(a:o)/^^(a;o) < en{xo)n^{xo) < 4>n > (98) 

which proves the first inequahty. 

Now, since 7^ < 2, 5^{1 - 5^)C^ < 6^{1 - 6^) < 1/4 and (ln8)/2 > 1, from the 
definition of we get that > 1 and prove the second inequahty. 

The last inequality is equivalent to 

and it follows from - S^) < 5^ < C^. 



4.2. Proof of Theorem [2771 We write 

^A-TXxn] = T^Ti{j^)^u[Tx\n I Tx\n < T^] + (1 - vr7^(^^))E^[T;v\7^ I '^A\7^ > , 
and notice that trivially 

0<K,[Tx^n\rx\n<T^]<T^- 



Moreover, by Theorem |2.4[ one gets 



T^ + :^{^-S*n}<^Arx\n\rxKii>T^]<T^ + :^{^ + S*n} 



Altogether this proves statement (i) . 
In the same way, we can write 



+ (1 - 7m{u))F,iT^T^ >^\r^^> T^) , 
If t > <PtiT^, the first term in the r.h.s cancel and we get 

P. (rA^7^ > 4) = (1 - v^7^(zv))P. (r;t^7^ > 4 I r;t\7^ > 



By Theorem 12 .41 we also have 



■rx\n > ^ I rx\Ti >Tt^) -e 



±_ 



that, together with the previous equality, provides the relation (ii) . 

4.3. Proof of Lemma [278l Let X denote the discrete version of X, such that X follows X 
at each ring of a Poissonian clock of intensity 1, and let N{t) denote the number of rings 
up to time t. By the natural coupling between the dynamics Xn and its discrete version 
X-ji we have 

k>l 

<T.I1 ^f^niXiiik - 1) = x)e^^{x)F{N{T^) > k) (99) 

k>l xeTl 
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5. Working with (k, A) -capacities 

5.1. Proof of the upper bound in Theorem 12.111 Let X denote the continuous time 
Markov chain on X defined, for k > 0,by the generator 

( k{f{x) - f{i)) iix = xe'iZ 

= <^ + K{f{x) - fix)) if X = X G 7^ . (100) 

[ iix = x£X\n 

This is a reversible process with respect to a measure jl defined as 

'^(^^ = I f ^(x) ifx = X e ^ • ^^^^^ 

Note that /i is not a probabihty measure. Let us denote by the harmonic measure on iZ 
associated with X\TZ, i.e., the probabihty measure on TZ defined by 

and with 

r K(x) = P^K < T;t^7e) ifx = xe7^ 

V^(x) = <1 ifx = xG7^ 

1^ if X = X G x\n 

With obvious notation, we then have 

E.,lf^] = ■ (103) 

Such kind of formula was introduced for the study of metastability in I116II . We refer to 
lecture notes ll3T]| for a derivation. 

Now setting u{x) = z>^(x) for all x eTZ, we can write 

and 

^ /x(K) + m ^ /i(K) + if^m 
c^{n,x\n) c^{n,x\n) c^{n,x\n) ' 

Inserting the above equalities in (I103D and multiplying by k, we then get 

Note that fi{TZ) , ij.{Vk) , 0^(11, X\TZ) and E;,[r;e\R.] do not depend on k. Then, in the limit 
of a vanishing k, it holds 

By Theorem |2.7[ this implies 

that is equivalent to the upper bound in (l35l) . 

For future use, we note that dividing Eq. (I104D by k, and then sending k to +oo, we get 
Together with (I105D . this implies ^ ^^x^) ~ c {n%R) ~ k equivalently 
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Lemma 5.1. 

c^{n,x\n) 

5.2. Proof of the upper bound in Theorem 12. 12i For any / e ^^(/i), we have 
Var^(/) = ^(Var^(/|l7e)) + Var^(//(/|l7e)) 

= ^(7e)Var^^(/|^) + ^(A'\7^)Var^^^(/|^^) ^-^^g^ 

+ //(7^)/.(A'\7^) (^7e(/|^) - /^A^7^(/^^)' ' 
Now, using the test function 



/ 



in the definition (l30l) of (k, A)-capacity, we get 

C^{n,X\n) < vCf) + K/.(7^)Var^^(4) + A^(A'\7^)Var^^^(/|^^) 

= {^^nU\J - t^mifi^jyX^if) + ^Km^^.M + V(^\^)Var^^,(/|^^)) , 

which provides an upper bound over (^fi'ji{f^^) — fJ.x\jz{f\x\Tz^^ " ■^PP^yi'^S that bound in 
Eq. (I108D . and from the definition oi (j)^{A, B),we get 
Var^(/) < /.(7^)Var^^(/|^)+/.(Af\7^)Var^^^(/|^^) 

+ (V{f) + K/.(7^)Var^^(/|^) + A/.(A'\7^)Var^^^(/|^^)) 0^(7^, A'\7^)"' (109) 



< -x#Q^U + max 



where in the last step we used that 

V{f) < fi{n)Vn{f\J + ^(A'\7^)P^7^(/|^^) • 
The upper bound in (|37l) follows directly. 

5.3. Proof of the lower bound of Theorem 12. Hi From inequality (1 109D applied to / = 
h*^ and A = +oo, and since = 0, we get 

Y^T^{h\) < M(7^)Var^^(/i^) + ^^^^||L^(7^)(l _ ^(7^)) {l + . (110) 

On the other and, by (fTOSl) . 

Var^(/i^) = A^(7^)Var^^(/l^) + /x(7^)(l - ^(7^)) . 

Inserting this formula in (IllOD . the term ^(7^)Var^^(/i^) cancels and then, dividing by 
/i(7^), we have 

'-''<^'^ cM7;) <'-'''^»{' + 4}- 

Now, modifying the transition probabilities along the edges that connect X\TZ toTZ with- 
out changing the (non-diagonal) ratios but modifying also /x, in order to keep a 
probability measure on X, we can send ^{IZ) to without changing fi-ji, and the ratio 
V{h^)/Ci^(lZ,X\R.). Then we obtain, for the modified system as well as for the original 
system, 

c^^{n, x\n) {1 + 4}"' < v{h*T^) . (112) 
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Now, dividing by /i(7^)||/i^||^ and using that, by Prop. 12.11 = Var^^(/i^) + 1 < 

1/(1 - e^), we get 



|l + ^r < „ =0^, (113) 



m(^) + " ^ ^'^i - ||/i^|||/x(7^) 

where the last equahty comes from Lemma [3TT] and the fact that 

Finally, using the convexity of the function x we obtain 

(1 - - ^1 ^ , (114) 

which concludes the proof of the lower bound in (l35l) . 

5.4. Proof of the lower bound in Theorem l2.12i We use the test function V^, for which 
we know that (see (l30l) ) 



^ \x\R. ^ 



so that 

V{V^) < C^{TZ,X\n) . (115) 
We then look for a lower bound on Var^(V^^). From (I108D we have 

var^(i/^) > ^(7^)^(A'\7^) (/i7e(KV) " ^^7^(v;\^)) ' ' 

thus we need to estimate f^n{V^\J and fix\n{VK\;^J- 
By the monotonicity in A, for all x G 7^ we get 

which implies, together with Lemma [5TT1 

^^^^{v^, i^J ^ -L ^M^)— ■ 

In the same way we have 
Altogether, we finally get 

6. Working with soft measures 

6.1. Proof of Lemma [2.13i If A = the first statement holds trivially since, in that case, 
^Tzx~^~ ^7?. a(^7?-,a)- If a > 0, we can write 

^^^'^Jr;^n,x <t) = Y.^^.'^JNn{t) > k){l - rn,xf'^^HiA^n,x) , 

k>l 

where N{t) is the number of clock rings inside TZ for the Poissonian clock associated to X. 
Taking the limit as t — ;> oo in the above equation, we get that 

1 = M7^,A(e7^,A)/0^,A (l^n,x = /^7^,A(e7^,A) , 

which provides identity i) . 

Let us now define the operator £^ on i'^ifj.n) as 

{c*T^^j){x) = -fix) + Pkx(^, y)f{y) G 7^ , / G eiim) (ii7) 
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and notice that, for any probability measure v on X, it holds 



/(XoG7e(t)l|._,>,|) 



vU^T^^^n. (118) 



The exponential law given in ii) follows from the above identity applied to v = and 

/ = 

Finally, since 1 — (/)^ ^ is a simple eigenvalue equal to the spectral radius of ^j^, for any 
x,y gTZ and in the large t regime, we have 

P,.(X o Gn{t = y) , Txyn,x > ~ c.,ii\^^{y)e-"^*^-^ , (119) 

where c^/U^ ^ is the canonical projection of 5x on the one-dimensional eigenspace associ- 
ated with Note that Cx is positive, otherwise we would have P^^ ^ (t;^\^,a > < 

(1 - //^ ^{x))e~^'^'^->^ + o(e~*'^^'^) while /^^ a > ^ri'i ^^is would contradict the exponential 
law. From ( I119D , and taking the limit of t — ;> oo, it follows 

lim ¥x{X oGn{t = y)\ Txsji,x > t) = /x^,;, . 

6.2. Proof of Lemma [2. 141 The result is once again a consequence of Perron-Frobenius 
theorem. Let xxiu) denote the characteristic polynomial of C*^x> ^^^^ '^^^ written as 
Xxiy) = {y- (p*Ti^x)0'iy)- If = iy- 4>*n^x)lhj) + ^(y) is the Euclidian division of a{y) by 
(y — <t'Tix)> then r{y) = r(0) / 0, and we have the Bezout identity 

- ^Qiy){y - 4>*n,x) = i • (120) 

In particular, for any x € TZ, :^6xa{C*^ x) = c^/x^ a the canonical projection of 6x on 
the eigenspace associated to j^i^ X' since Cx > as previously noticed, we have 

^'R.,x = ^ ^~TP — u~ ■ '^^^^■^ 

Since a{y) = (J^^*^ y the above equation expresses the map A i-^ a ^ composition 
of continuous functions of A. 

6.3. Proof of Proposition [27151 As far as 4'n,x is concerned, continuity and monotonicity 
derive from continuity and monotonicity of eT^^AC^;) for any x G TZ. We then consider the 
other parameters. The continuity follows from the continuity of the eigenvalues as root of 
the characteristic polynomial. To prove the monotonicity, we notice that when A decreases 
to zero, p^x{^,y) grows for for all x and y in 7^ as well as c-ji^\{x,y) for any distinct 
x,y gTZ. From the variational characterization of 0^ x, i.e. 



'-"R.,x 



mm 



{ (/, -^n,xf)n : (/, /)7^ = 1 , / > o} (122) 



= min_ Yl fix) - Y.PnA^^y)f(y) (123) 

'^f'fj^^-^ x,y<^n \ yen J 

where the restriction / > comes from the fact that, by the Perron-Frobenius theorem, the 
right eigenvector has positive coordinates, we see that cp"^ ^ is decreasing in A. Similarly 
using 



77e,A = mm < 



\ Y cn,x{x,y){f{x) - f{y)f : Var^^(/) = 1 ^ , (124) 
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we see that 77^ ^ is increasing in A. As a consequence ^ is decreasing in A, and we have 

^■R,o = — - = — ■ = . (125) 

77^,o m,o 



6.4. Proof of Theorem [27T71 Proof of ([55]): We first write 



,{X{ts) = x\X{ts) GTZ) 



1 



\{X{Ts) G 7^) 



^Y.^-('o>i,Xin) = Xi) (126) 



i>0 XiGX 



Now, conditioning on a^, and setting PJ'^ = 1 ^k), we get 



P,(X(r5) =x|X(r5) e7^) 



1 



- AX{Ts)en) ^^^^^ 



X E 



'i^(Ts)en) ^^^^^ 



(127) 



where the second equality comes from the independence between X, and ax. Since 



i>0 XiSA' 

from (I127D we get 



(128) 



.(X(r5) = X I X(t5) e 7^) 



1 



1 



X E 



^ ^ F,{io>i,XiTi[ 



Xi 



i>0 Xi&X 

o Gni'yn) =x\a^< r;t\7^,A) 
^7e,A 



(129) 



An analogous expression can be found for ^•'^-^'^'^^^ x\x{ts)&x^) _ -j^ r^^^ result then follows 

from Th J2.16[ and in particular for the equivalent of Th. 12. 4[ 
To prove inequality (l56l) . we first state the following lemma. 

Lemma 6.1. Let T > and {(j, : i > 1} be a sequence of independent exponential random 
variables of rate k such that e'^^ — 1 < 1. if = min{i > 1 : fj, > T}, then 



(130) 
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Proof of Lemma \6l] Using the property of the exponential distribution, we have 

> ^ ) = J^P(A^ = n)P [ ^(Ji > ^\ai<T,...,an-i<T,an>T 

=1 '^J n>l \j=l J 

< ^¥{N = n)P ( cJn > - - (n - 1)T \an > Tj 

n>l \ ^ / 

= ^ P(A^ = n)P ( > - - nT ) (131) 



n>l 



Ed 



g-KT\n-lg-KTg-t+nKr 



ra>l 

-t 



e 



n>l ^ ^ 

which concludes the proof. □ 

Coming back to the proof of Th. |2.17[ we first notice that if ts > + j), then L'ji{ts) > 
i or Lx\jz{ts) > ^. As a consequence, defining 

A-jz = {KL-jiirs) V XLxynin) = kLti{ts) > t} 
Ax\n = {kLtz{ts) V XLx^nirs) = ^Lxyni^s) > t] 

so that ^{A-ji) + P(^A\R.) < 1, we have 

P.(r5 > t(i + {)) = P.(t5 > t(i + {)\ An)^MTl) + ^uirs > + i)| A;t^7^)IP.(^A^7^) • 

Using the independence between a^,, a\ and X, together with the previous lemma, we 
finally get 

P, {T-6>t{^^ + ^))<e-'{j^^]{K{An)+K{Ax^Tz))<e-'{j^^] . 

6.5. Proof of Theorem l2.20l Considering the extended electrical network associated with 
C^(^, B), the proof of the upper bound is similar to that of Theorem 12.111 and we omit it. 

As far as the lower bound is concerned, the situation is much more delicate. In the proof 
of the analog Theorem [2TTT] we used a partial Poincare inequality to control the mean exit 
time from IZ. Here we will do the same, with the difference that we will have to work 
on the whole space X and not only on IZ. Since ^ is concentrated on IZ, we will first 
compare it with another quasi-stationary measure, jl*^, that spreads on the whole space 
X. Then we will control the escape rate of /i^ from X with the spectral gap estimated in 
Theorem 12. 12[ This second step will be the easier one. 

Let fi*^ be the quasi-stationary measure on X associated with the Markovian process X 
on X = X U xyiZ with generator C defined, for some A > 0, by 

r ifx = xE7^ 

{Cf){i) = \ {Cf){x) + X{f{x)- fix)) iii = xeX\R . (133) 
[ A(/(x) - f{x)) if X = X G X\n 

With obvious notation, for any x ^ X and any probability measure v on X, we have 

A:^(x)=limP,(X(t) = x|f >f). (134) 

We want to show that Jj*^ is in some sense close to ^j^. The difficulty comes from the fact 
that while Jl*^ is defined through the global time of the dynamics, /z^ ^ is defined through 
the local time in IZ. In particular, it is not true that //^ ^ = fi*p^{- \ TZ) . 
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For any a > 1 and T > 0, there is a large enough t > T such that, for all x 

^F,,^jX{t) = x,f—>t) 



a 



P.* (f > t) 

Recalling the definition of the times Tj given at page [T2l we set 

r = min{ri > : Tj > t - T} . 
The numerator in (I135D then satisfies 

P,^ = X, f_ > t) > P,^,,(X(t) = X, f_ > t > T) 

^/^^,A {^{x(t)=x,f—>t}\^r\ l{Mf^>r} 



(135) 



X\TZ 



E 



'xiT){X{t -T) = x,f^>t- r)l{Mf_>r} 



To what concerns the denominator, we have 

P.. if >t)<p.. if >tAr) = Pu* (r<tAf )+p.. {T^f >t).(i36) 



Now, since P^^ JX(ti) = • | f_ > n) = /i^^^, we get 



< 



Now, similarly to the proof of Theorem (I2.17D but, since we do not care about condi- 
tions (I53H54D . without the correcting factor we have, on one hand. 



Vi > T) < e"^('^ + 9 



(137) 



and, on the other hand. 



\*^^{t, < T, X{n) G 7^) = f,*^ JX{n) G 7^) - F^^^^in > T, X{n) e U) 



= 1 

> 1 

> 1 



(138) 



where the last inequality comes from the following direct computation for k / A and 



+ ex. > T) = < = e'^^^^^ 



and a continuity argument for the case k = X. 
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We conclude that 



< (T <tA f ] 



1 + - 
V 1 



(139) 



< 



(r < t A f- 



and coming back to (fT35]) . with ai = l- - e ^V^^a 



1 , 1 



ai 



X 



irmt - r) = X, f_ > t - r)l|Mf_>r} 



p„. (r < t A f. 



(140) 



Now we note that starting from jjL^ ^, and conditionally on {T < tA t— }> the law of X{T) 
is close to fi^ in total variation. Indeed, 

= P^;^ jX(ri) = • I x(ti G 7^), n < T)F^*^Jn <T\x{ri G 7e)) (141) 

+ w^^^ jxin) = ■ I x{n en),n> t)¥^*^ Jti > T \ x{ti g 7^)) 

and, by (fT38]) 

-;.,A -rr '1^-^ 



P^;^,A(^l < r| ^(n) G 7^) > P^^ < T, X(ri) G 7^) > 1 
As a consequence 



— e 



^.;,,,(^(r) = -|r<tAf_)-^^ 



< 



1 — Ql 



TV ai 



(142) 



This is useful because, from the estimate 

Fx(T){Xit -T) = x,f—>t-T)> FxiT){X{t - T) = X, r;t\7e,A >t-T) 

> Px(r)(^(i -T) = x\ Tx\n,x > t - r)Px(r)(^A\7^,A > T) 



and since 



{tx\tz,x >T) = e '^T^'^ , we get that for all a > 1, 

1 — Ql ' 



Ql 



(143) 



u{x) = ^v{x) , 

where v{x) is a probability measure such that — ;^||tv < ^^f^- particular, 

vr7e,A(z^) < vr7e,A(^^,A) + ^ = 1 - e-'^-A^i.A + 1^ = i _ ^3 , 
where, we recall, TXn.xi'^^) = P^/(''"A\7^.A < T^), and by Theorem ( I2.19D we get 

lE^tW = ^^[^:^ + '^a] ^ #7 {1 + '^^-A^^^.A - -^^.a} + i = ^ (144) 

Finally, dropping the coefficient 1/a in (I143D since the inequality holds for all a > 1, we 
get 



(145) 



n,x 
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This gives a lower bound on (p^ ^, provided we have an estimate on c^^. We get such an 
estimate by comparison with the spectral gap. By Lemma 13.11 applied with TZ = X and, 
with obvious notation, f = h\, we have 



> 



XI 



\h* l|2 



1 V{h 



X) 



\h* l|2 



|/i^||2 Var^(/i^) 



> 1 



1 



\h* l|2 



7- 



Now, 



^x\\ 



> 



f^x( 



fl[X) 



(.«l«2j' 

/x(7l) 



so that, by Theorem [ZT2l 



''X 



> 1 



{aia2fj /i(7^)(l-M(7^)) i + niaxf^,^' 



(146) 

(147) 
(148) 

(149) 
(150) 

(151) 



Plugging this last estimate in (I145D and using the convexity of x 1-^ 1/(1 + x), we get the 
desired result. 



6.6. Proof of Theorem [2:211 By Theorem [27T61 for all x in X, 



max 



•x{x{T) = y 



- '^A^7^,0• 



Then 



1 

2 
1 

2'^x<ji, 
1 



(152) 

(153) 
(154) 
(155) 



Also ||/U;v\-R.— Aill = m(^) and the upper on ||z^a;— follows from the triangular inequality. 
Now, for all t > 0, 



Wx{x{t) = ■)- /iiixv < hx - /uIItv + Px-(r > t) 

and to prove our mixing time estimate, it is sufficient to show 



\\T> 



1 + 7/ ^ 1 ^ 



(156) 



(157) 



To obtain such an estimate we give an upper bound on the mean value of T and use Markov 
inequality. 

With T' = ax A Til, conditionally to X(7'') £U,we have 



(158) 
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Then 



(159) 
(160) 



and 



r > 



(1 - 



< ^ — + ^^ + ^^['Pn,xTTi,x + OTi,\) (161) 



< ^ + o I ^ ^ <P-77 i-iT? A + 0% 



2 

1 + ?? 



2 V A 



(162) 
(163) 



Appendix A. A very crude estimate 



We prove here lemma [231 The first inequality is obvious from the definition of Q^. Let 
X denote the discrete version of X, such that X follows X at each ring of a Poissonian 
clock of intensity 1, and let N{t) denote the number of rings up to time t. Then, for z eTZ, 
we have 

= Urn P,.(X(t) = z I r;t^7^ > 



= lim V F^iXik) = z\fx\n> k)¥{N{t) = k) 

fc>0 

> lim ^o.iX{k + Dn) = z\fx^>k + L»7e)F(iV(t) = + Z^Te) (i64) 



fc>0 



fc>0y67^ 



X P(iV(i) = k + Dn) 



where we used the notation for the hitting time of the chain X on X\1Z. Since for all 
y eTZwe have 



,{X{Dn) = z\fx\n> Dn) > ^y{XiDn) = z,fxyn> Dn) > e"^'^^^ 



we get 



f-iniz) > e-^^^'^ lim E ^-(^(^) =y\^m> k)nN{t) = k + Dn) 

fc>0ye7^ (165) 



g-zi7ji^7j lini P(7V(t) > Dn) = e.-'^T^'^^ . 



References 

[1] W. Thomson, P. G. Tait, Treatise on Natural Philosophy, Oxford (1867). 

[2] J. N. Darroch, E. Seneta, On quasi-stationary distributions in absorbing discrete-time finite Markov 

chains, J. Appl. Probability 2, 88-100 (1965). 
[3] J. L. Lebowitz, O. Penrose, Rigorous treatment of metastable states in tlie van der Waals-Maxwell 

Theory, J. Stat. Phys. 3, 211-241 (1971). 
[4] A. D. Wentzell, Formulas for eigenfunctions and eigenmeasures that are connected with a Markov 

process, Teor. Verojatnost. i Primenen. 18, 329 (1973). 
[5] M. Cassandro, a. Galves, E. Olivieri and M. E. Vares, Metastable behaviour of stochastic dynamics: 

a pathwise approach, J. Stat. Phys 35, 603-634 (1984). 
[6] M. I. Freidlin, a. D. Wentzell, Random Perturbations of Dynamical Systems, Springer-Verlag (1984). 



METASTABILITY AND QUASI-STATIONARY MEASURES 



31 



[7] E. J. Neves, R. Schonmann, Critical droplets and metastability for a Glauber dynamics at very low 

temperature. Comm. Math. Phys. 137, 209-230 (1991) 
[8] E. J. Neves, R. Schonmann, Behaviour of droplets for a class of Glauber dynamics at very low temper- 
ature, Prob. Theory Relat. Fields 91, 331-354 (1992). 
[9] E. SCOPPOLA, Metastability for Markov chains: a general procedure based on renormalization group 
ideas, in ProbabiUty theory of spatial disorder and phase transition Cambridge (1993) Ed. G.R.Grimmett 
- Kluwer Acad.Publ.(1994) 303-322. 
[10] E. ScoppOLA, Renormalization and graph methods for Markov chains, in Advances in Dynamical Systems 

and Quantum Physics S.Albeverio, R.Figari, E.Orlandi, A.Teta Ed. - World Scientific 1995. 
[11] G. Ben Arous, R. Cere, Metastability of the three-dimensional Ising model on a torus at very low 

temperature. Electron. J. Prob. 1 (1996). 
[12] P. Dehghanpour, R. Schonmann, Metropolis dynamics relaxation via nucleation and growth. Comm. 

Math. Phys. 188, 89-119 (1997). 
[13] J. R. NORRis, Markov Chains, Cambridge University Press (1997). 

[14] P. Mathieu, p. Picco, Metastability and convergence to equilibrium for the random field Curie-Weiss 

model, J. Statist. Phys. 91, 679-732 (1998). 
[15] R. Schonmann, S. Shlosman, Wulff droplets and the metastable relaxation of kinetic Ising models. 

Comm. Math. Phys. 194, 389-462 (1998). 
[16] A. BoviER, M. ECKHOFF, V. Gayrard AND M. KLEIN, Metastability and low lying spectra in reversible 

Markov chains, Commun. Math. Phys. 228, 219-255 (2002). 
[17] A. BoviER, F. Manzo, Metastability in Glauber dynamics in the low temperature limit; Beyond expo- 
nential asymptotics J. Statist. Phys. 107, 757-779 (2002). 
[18] G. Ben Arous, A. Bovier and V. Gayrard, Glauber dynamics of the random energy model. 1. 

Metastable motion on the extreme states, Commun. Math. Phys. 235, 379-425 (2003) 
[19] A. Bovier, M. Eckhoff, V. Gayrard and M. Klein, Metastability in reversible diffusion processes 1. 

Sharp estimates for capacities and exit times. J. Eur. Math. Soc. 6, 399-424 (2004). 
[20] W. HuisiNGA, S. Meyn, C. Schutte, Phase transitions and metastability in Markovian and molecular 

systems, Ann. Appl. Prob. 14, 419-458 (2004). 
[21] M. Jerrum, J. B. Son, P. Tetali and E. Vigoda, Elementary bounds on Poincare and log-Sobolev 

constants for decomposable Markov chains, Ann. Appl. Prob. 14, 1741-1765 (2004). 
[22] A. Bovier, V. Gayrard and M. Klein, Metastability in reversible diffusion processes. 2. Precise esti- 
mates for small eigenvalues J. Eur Math. Soc. 7, 69-99 (2005). 
[23] A. Gaudilliere, E. Olivieri and E. Scoppola, Nucleation pattern at low temperature for local Kawasaki 

dynamics in two dimensions, Markov Processes Relat. Fields 11, 553-628 (2005). 
[24] E. Olivieri, M. E. Vares, Large deviations and metastability Cambridge University y Press (2005). 
[25] A. Bovier, F. den Hollander and F. Nardi, Sharp asymptotics for Kawasaki dynamics on a finite box 

with open boundary conditions Probab. Theor Rel. Fields. 135, 265-310 (2006). 
[26] A. lOVANELLA, B. ScoppOLA AND E. ScoppOLA, Some spin glass ideas apphed to the chque problem, J. 

Stat. Phys. 126, 895-915 (2007). 
[27] O. Bertoncini, J. Barrera M. and R. Fernandez, Cut-off and exit from metastability: two sides of 

the same coin, C. R. Math. 346, 691-696 (2008). 
[28] A. Bovier, A. Faggionato, Spectral analysis of Sinai's walk for small eigenvalues Ann. Probab. 36, 

198-254 (2008). 

[29] E. N. M. CiRiLLO, F. R. Nardi AND C. Spitoni, Metastability for reversible probabihstic cellular au- 
tomata with self-interaction, J. Sfat. Phys. 132, 431-471 (2008). 

[30] A. BiANCHi, A. Bovier and D. Ioffe, Sharp asymptotics for metastability in the Random Field Curie- 
Weiss model, EJP 14, 1541-1603 (2009). 

[31] A. Gaudilliere, Condenser physics applied to Markov chains - A brief introduction to potential theory, 
arXiv:0901.3053 

[32] A. Gaudilliere, F. den Hollander, F. R. Nardi, E. Olivieri and E. Scoppola, Ideal gas approxima- 
tion for a two-dimensional rarefied gas under Kawasaki dynamics, Stoch. Proc. and Appl. 119 737774 
(2009). 

[33] J. Beltran, C. Landim, Tunneling and metastability of continuous time Markov chains J. Stat. Phys. 
140, 10651114 (2010). 

[34] L. MiCLO, On absorption times and Dirichlet eigenvalues, ESJAM ; PS 14, 117-150 (2010) 

[35] A. Bovier, F. den Hollander and C. Spitoni, Homogeneous nucleation for Glauber and Kawasaki 

dynamics in large volumes and low temperature Ann. Prob. 38, 661-713 (2010). 
[36] A. Bianchi, A. Bovier and D. Ioffe, Point-wise estimates and exponential laws in metastable systems 

via coupling methods, Ann. Prob., in press, available at arXiv:0909.1242 
[37] J. Beltran, C. Landim, Metastability of reversible condensed zero range processes on a finite set, 

|arXiv:0910.4089l 



METASTABILITY AND QUASI-STATIONARY MEASURES 



32 



[38] J. Beltran, C. Landim, Metastability of reversible finite state Markov processes, arXiv: 1009.40901 
[39] A. Gaudilliere, B. Scoppola, E. Scoppola and M. Viale, Phase transitions for the cavity approach to 

the chque problem on random graphs, arXiv:1011.2945 
[40] R. Cere, F. Manzo, Nucleation and growth for the Ising model in d dimensions at very low temperatures, 

arXiv:1102.174 

A. BiANCHi, Dip. di Matematica - Universita di Bologna, Piazza di Porta S. Donato 5, 40126 
Bologna, Italy 

E-mail address: alessandra.bianchi7@unibo . it 

A. Gaudilliere, LATP / IMPA, Universite de Provence, CNRS, 39 rue F. Joliot Curie, 13013 Mar- 
seille, France 

E-mail address: gaudilli@cmi.univ-mrs.fr 



